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The degeneracy between star-formation parameters in 
dwarf galaxy simulations and the Mgtar — ^haio relation 
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CN ■ 1 INTRODUCTION 

J!* Dwarf galaxies are the most common type of galaxy in the 
• 1^ " local universe but also the faintest and least easy to observe. 
^^ ] In the ACDM cosmology, our universe consists of matter, 
^ n both luminous and dark, and dark energy, which is respon- 
. . .' sible for the accelerating expansion of the universe. Galaxies 
form when gas collapses in dark matter halos. Baryons, be 
it in the form of gas, dust or stars, are the most accessible 
form of matter, emitting radiation over the whole electro- 
magnetic spectrum. Dark matter, on the other hand, as it 
only interacts gravitationally, is much more difficult to "ob- 
serve" . 

There have been many attempts to estimate dark halo 
masses and mass-to-light ratios for galaxies and clusters of 
galaxies from direct observations. The se include methods 
that make use of gravitationa l lensing (JMandelbaum et al.l 
bood : iLiesenborgs et al.l |2009|) , dynamical modeling of the 



ABSTRACT 

We present results based on a set of N-Body/SPH simulations of isolated dwarf galax- 
ies. The simulations take into account star formation, stellar feedback, radiative cool- 
ing and metal enrichment. The dark matter halo initially has a cusped profile, but, at 
least in these simulations, starting from idealised, spherically symmetric initial con- 
ditions, a natural conversion to a core is observed due to gas dynamics and stellar 
feedback. 

A degeneracy between the efficiency with which the interstellar medium absorbs 
energy feedback from supernovae and stellar winds on the one hand, and the density 
threshold for star formation on the other, is found. We performed a parameter survey 
to determine, with the aid of the observed kinematic and photometric scaling relations, 
which combinations of these two parameters produce simulated galaxies that are in 
agreement with the observations. 

With the implemented physics we are unable to reproduce the rela tion between 
the stellar mass and the halo mass as determined by iGuo et ahl (|2Q1Q[ ) , however we 
do reproduce the slope of this relation. 
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observed propertie s of a kinematical tracer such as stars or 
plane t ary nebulae (IKronawitter e t al. 2QQ d:lDe Riicke et al.l 
I2OO6I : iNapolitano et al.1 l201l'i nBarnabe et al.1 l2009l ). One 
thing virtually all these works have in common is the rel- 
atively limited s ize of the data set they are based on. 
IGuo et al.l (|20ld l determined the halo mass as a function 
of stellar mass for a large sample of galaxies using a statis- 
tical analysis of the Sloan Digital Sky Survey, which yields 
the stellar masses, and the Millennium Simulations, which 
yield the dark-matter masses. In the range of the most 
massive halos and bright galaxies, the derived Mstar-Mhaio 
relation, which is of the form Mstar ex: M^afo? is found 
to be in good agreement with gravitational lensing data 
([Mandelbaum et al.ll2006l V Below a halo mass of Mhaio ^ 
10^^'^ M(T), this r elation becomes much steeper: Mstar oc 
M]^'^^. IGuo et al.l (|2010) extrapolate the latter relation into 
the dwarf regime, where Mhaio ^ 10^° Mq. This leads 
then to the prediction that faint dwarf galaxies with stel- 
lar masses of the order of Mstar '^ 10^ Mq should live in 
comparatively massive Mhaio ^ 10^° Mq dark-matter halos. 

The IGuo et al.1 (l20ld ) Mstar-Mhaio relation was 
compared with that f ound in simulations of dwarf 
galax ies (IValcke et al.l 
2OO9I: iGovernato et al 
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ISaw ala et alj (|201lh . They found that simulated dwarf 
galaxies had stellar masses that were at least an order 
of magnitude higher at a given halo mass than predicted 
by IGuo et alJ (|20ld V There could be several causes for 
numerical dwarf galaxies to be overly prolific star formers: 



• The star formation efficiency could be too high be- 
cause of an under estimation of the feedback efficiency. 
IStinson et alJ (|20Q6h investigated the influence of the feed- 
back efficiency on the mean star formation rate (SFR). The 
general trend they have observed was a decrease of the mean 
SFR when increasing the feedback efficiency. 

• IStinson et al.l (|2006|) also reported finding a decreasing 
mean SFR with increasing density threshold for star for- 
mation. Recently, high densi ty thresholds for s tar fo rmation 
have come in vogue, see e.g. iGovernato et al.l (|20ld V 

• Dwarf galaxies, due to their low masses, are expected to 
be particularly sensitive to reionisation. Not properly tak- 
ing into account the effects of reionisation may lead to an 
overestimation of the gas content of dwarfs and an underes- 
timation of the gas cooling time. 

• Dwarf galaxies are metal poor and hence also dust poor. 
This lowers the production of H2 niolecules and caus es poor 
self-shielding of molecular clouds (jBuvle et al.ll2QQ6l ) which 
could be expected to inhibit star formation. Not taking these 
effects into account will lead to an overestimation of the SFR 
(|Gnedin et al.ll2QQ9l ). 



Using the high values for the density threshold above 
which gas particles become eli gible for star formation, de- 
noted by psF, as promoted bv IGovernato et al.l (|2010|), in 
combination with radiative cooling curves that allow the gas 
to cool below 10^ K (Maio et al.ii2QQ7i ). makes the gas col- 
lapse into small, very dense and cool clouds before star for- 
mation ignites. If the supernova feedback epB, defined as the 
fraction of the average energy output of a supernova that is 
actually absorbed by the interstellar medium (ISM), is too 
weak to sufficiently heat and/or disrupt such a star-forming 
cloud, one can consequently expect the mean SFR to be very 
high, leading to overly massive (in terms of Mstar) dwarfs. 
Therefore, one could hope to remedy this situation by in- 
creasing cfb accordingly. In that case, a correlation between 
epB and psF would be expected to exist. 

In the present paper, we analyze a large suite of nu- 
merical simulations of isolated, spherically symmetric dwarf 
galaxies in which we varied both the feedback efficiency cfb 
and the density threshold psF- Our goal is to investigate (i) 
if such a correlation between cfb and psF exists and, if it ex- 
ists, how to break it, (ii) which epB/yOSF-combinations lead 
to viable dwarf galaxy models in terms of the observed pho- 
tometric and kinematic scaling relations, and (in) how well 
these models approximate the aforementioned Mstar -Mhaio 
relation. 

In section [21 we give more details about the numerical 
methods that are used in our code. An analysis of the simu- 
lations is given in section [3l where some details are given of 
the NFW halo that is used for the simulations and a large 
set of scaling relations are plotted comparing our models to 
observations. In section [4l we discuss the obtained results 
and conclude. 



2 NUMERICAL DETAILS 

We use a modified version of the Nbody-SPH code Gadget- 
2 (Springel 200^)- The original Gadget-2 code was ex- 
tended Math_£ta^formation, feedback and radiative cool- 
ing by IValcke et al.l (|2008l ) . While the initial conditions of 
the simulations are cosmologically motivated (see below) , we 
do not perform full cosmological simulations. Our approach 
yields a high mass resolution at comp aratively low compu - 
tational cost . Still, prev ious work by IValcke et al.l (|2008l l, 
IValcke et al! (|2Q1Q| ) and ISchroven et al.1 (|2Qllh lias shown 
that with this code realistic dwarf galaxies, following the 
known photometric and kinematic scaling relations, can be 
produced. We set up the simulations using 200,000 gas par- 
ticles and 200,000 DM particles. Depending on the model's 
total mass, this results in gas particle masses in the range 
of 350 — 2, 62OM0 and DM particle masses in the range of 
1, 650 — 12, 38OM0. We use a gravitational softening length 
of 0.03 kpc. 

Our results are visualized with our own software pack- 
age HYPLOT. This is freely available from SourceForg4j 
and is used for all the figures in this paper. 



2.1 Initial conditions 



Our models ar e set up, as in IValcke et al.l (|2008l . |201QI 1: 
ISchroven et al.l (|201lh , with a spherically symmetric dark 
matter halo and a homogeneous gas cloud. This gas cloud 
has a density of 5.55 pcrit, with pcrit the critical density of 
the universe at the halo's formation redshift, here taken to 
be Zc = 4.3. This is equivalent with a number density for 
the gas of 0.0011 hydrogen atoms per cubic centimeter. We 
use a flat A-dominated cold dark matter cosmology with the 
following cosmological parameters: h = 0.71, Qtot = 1, ^m = 
0.2383, Qdm = 0.1967. The baryonic mass fraction will be 
the difference between Qm and Qdm, in practice it will have 
a value that is 0.2115 times that of the dark-matter. At the 
start of the simulations the gas particles are initially at rest, 
their initial metallicities are set to 10~^ Zq and their initial 
temperature is 10^ K. The dark m atter halo has a NFW 
density proflle (JNavarro et al.|[l996l l: 



PNFw(r) 



(r/rs)(l + r/rs 



(1) 



where ps and rg are, respectively, the characteristic den- 
sity and the scale radius. In order to flx the values of 
these par ameters, we use the corre lation betwee n them 
found by IWechsler etHI (I2OO2I ) and iGentile et all (120041 ). 
which makes the NFW density distribution essentially a one- 
parameter family of the dark matter virial mass, Mdm- The 
relations we use for ps, Ts and the concentration parameter 
c (=rmax/rs) are : 
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^ http://sourceforge.net/projects/hyplot/ 
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Here, rmax is the halo's virial radius. At rmax, the DM halo is 
truncated and the density drops to zero, so the entire mass 
Mdm is situated inside the radius rmax- 



Table 1. Details of the basic spherical dwarf galaxy models that 
were used in the simulations. Initial masses for the DM halo and 
gas are given in units of lO^M©, radii in kpc. 



model Mdm,! Mg 



2.2 Criteria for Star formation 

Star formation is assumed to take place in cold, dense, 
converging and gravitationally unstable molecular clouds 
(JKatz et al.l lll996). Gas particles that fulfill the star forma- 
tion criteria (SFC) are eligible to be turned into stars. These 
SFC are: 



Pg ^ PSF 

T ^ Tc = 15000K 
V.v ^ 0, 



(5) 
(6) 
(7) 



with pg the gas density, T its temperature and v its velocity 
field. psF is the densi ty thre shold for star formation. We em- 
ploy a Schmidt law ([Schm idt 1959. ) to convert gas particles 
that fulfill the SFC into stars: 



dps 
dt 



dt 



c ^ 



(8) 



with ps the stellar density and c^ the dimensionless star for- 
mation efficiency. The timescale tg is taken to be the dynam- 
ical time for the gas 1/ y/^nGpg. Here, we choose c^ = 0.25. 
IStinson et al.l ((20061) showed that the influence on the mean 
SFR of the value of c^ with values in the range of 0.05 to 
1 is negligible. Lowering c^ reduces the star formation effi- 
ciency as well as the amount of supernova feedback, causing 
more particles to fulfill the density and temperature crite- 
ria. This compensates for the lower value of c^ , producing a 
SFR which is roughly independent of c^. 

1 i "T^N 

evaz et al.l (|2009l ) also investigated the infiuence of c^ 
by varying it between the values of 0.01 and 0.3. They con- 
cluded that the star formation history is mainly determined 
by the initial total mass with a minor infiuence of c^. Self- 
regulating models, in which star formation occurs in recur- 
rent bursts due to the interplay between cooling and super- 
nova feedback, were achieved for c^ ^ 0.2. Such models best 
resemble real dwarf galaxies. 



2.3 Feedback 

We consider feedback from star particles by supernova la 
(SNIa), supernova II (SNII) and stellar winds (SW). They 
deliver energy and mass to the ISM and enrich the gas. Feed- 
back is distributed over the gas particles in the neighborhood 
of the star particle according to the SPH smoothing kernel. 
Each star particle represents a single-age, single-metallicity 
stellar population (SSP). The stars within each SSP are dis- 
tributed according to a Salpeter initial mass function : 



^{m)dm = Am 



(9) 



with X = 1.35 and A = 0.06. The limits for the stellar masses 
are mi = 0.01 M© and ttIu = 60 Mq. The energy release 
of a SN is set t o ^tot = 10^^ erg an d that by a SW to 
Etot = 10^° erg (JThorntonet al.lll998l ). The actual energy 
injected into the ISM is implemented as epB x ^tot, where 
epB is a free parameter. 



N03 
N05 
NOG 
NOT 

N08 
N09 



330 
660 

825 
1238 
1654 
2476 



70 
140 
175 
262 
349 
524 



0.412 
0.566 
0.627 
0.756 
0.863 
1.040 



17.319 
21.742 
23.393 
26.755 
29.428 
33.634 



2.4 Cooling 

Metallicity-dependent radiati ve cooling is implemented us - 
ing the cooling curves from [Sutherland Sz Dopital (| 19931 ). 
With this recipe it is possible to cool gas to a minimum 
tempe rature of 10^ K. We also implemented the lMaio et al.l 
(| 2 00 71 ) cooling curves, making it possible for particles to cool 
below 10^ K. 



2.5 Production runs 

In Table [T] we give an overview of the parameters that were 
used to set up the models. A benefit of our code is that we 
can retain the same initial conditions and easily adapt our 
parameters to perform a detailed parameter survey. In the 
remainder, we will quantify the density threshold by nsF 
expressed in hydrogen ions per cubic centimeter (so psF — 
1 amuxrisr). At the start of the simulations, the models only 
contain gas and dark matter. During the first few 10^ years, 
the gas collapses in the gravitational potential well of the 
DM. The simulations run for 12.22 Gyr, tifl z = 0. 

In the literature, a larg e variety of va l ues fo r the den- 
sity threshold can be found. IStinson et al. (|2006r) use a low 
density threshold of 0.1 cm~^ while iGovernato et al.l (|2010h 
use a high density threshold of 100 cm~^ which, these au- 
thors argue, is a better representation of the conditions in 
star- forming r egion s in real galaxies. The simulations of 
ISawala et al.l (|201lh have been performed with a density 
threshold of 10 cm~^. In this paper, we increase the den- 
sity threshold from nsF = 0.1 cm""^, over nsF = 6 cm~^ to 
^SF = 50 cm~^. For the fiducial series of low-density thresh- 
old simulations, we matched the nsF = 0.1 cm""^ with a feed- 
back efficiency of cfb =0.1. For the intermediate-density 
threshold simulations, with nsF = 6 cm~^, we varied the 
feedback efficiency between cfb =0.1 and 0.9. Finally, for 
the high-density threshold simulations, with nsF = 50 cm~^, 
we varied the feedback efficiency between cfb =0.3 and 0.9. 



3 ANALYSIS 

3.1 The NFW halo 

The DM halo is constructed using a Monte Carlo sampling 
technique. First, for each particle, the three position coor- 
dinates in spherical coordinates (r, ^, 0) are generated, r 
is drawn from the density profile pnfw using a standard 
accept ance-rejectance technique, cj) and cos(^) are drawn 
from uniform distributions over the intervals [0, 27r] and 
[—1, 1], respectively. Next, Vr^ ve and -u,^ are drawn from the 
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0.040 
0.035 
0.030 



Figure 1. The density profile of the N03 NFW halo for different 
simulations: in the upper panel only DM was included, in the 
central panel DM, gas is included but star formation was turned 
off. The bottom panel shows the results of a simulation with DM, 
gas and star formation. 



isotropic distribution function for the NFW model, again 
with an accept ance-rejectance technique. This isotropic dis- 
tribution function was constructed from the N FW density 
profil e using the standard Eddington formula (|Buvle et alj 
120071 1. For each particle, a symmetric partner was con- 
structed with position coordinates (r, — ^, —cj)) and velocity 
coordinates (— I'r, — I'e, — i'^). This drastically improved the 
stability of the central parts of the halos. The very inner 
part of the steep cusp of the NFW model is populated by 
relatively few particles, destroying its spherical symmetry 
and introducing unbalanced angular momenta. This initial 
deviation leads to the ejection of particles from the cusp 
and triggers a more widespread dynamical response of the 
DM halo, over time erasing the inner cusp. Introducing the 
partner particles, cancelling out the angular momenta and 
increasing the symmetry of the particles' spatial distribu- 
tion, greatly alleviates these problems. Such techniques for 
constructing "quiet" initial condit ions have been applied be- 
fore with great success, see e.g. ([Sellwood k, Athanassoulal 
1 19861 ). The improvement of the stability of the DM halo in 
simulations with a "quiet" start over simulations without a 
"quiet" start is illustrated in the top panel of Fig. [T] where 
the density distribution of both haloes at 2; = is plotted 
as red and green dots, respectively. 

First, to test the stability of the NFW halos, we ran 
several simulations for the N03 and N05 mass models: 

Run 1: only DM 

Run 2: DM and gas but no star formation 




0.000 
200 



nsF = 0.1 cm , epB = 0.1 
nsF = 6 cm~^, epB = 0.1 



nsF = 6 cm , epB = 0.3 
nsF = 50 cm"^, epB = 0.3 




4 6 

time (Gyr) 

Figure 2. Top panel: the SFR of several NOT models as a function 
of time. Bottom panel: the stellar mass as a function of time. 



Run 3: DM and gas and star formation 

For t hese test simula tions an nsF of Q-1 cm ~'^ (|Katz et al.l 
Il996l ) and cfb of 0.1 (JThornton et al.lll998l ) was used. 

Fig. [1] shows the density profile of the test simulations 
for the N03 mass model. From the upper panel, it is evident 
that the DM density of the DM-only simulation remains 
stable and cusped until the end of the simulation. The sim- 
ulations presented in the middle and bottom panels, show a 
clear conversion of the cusp into a core over time. Moreover, 
the width of the core depends on the mass of the system, 
with more massive halos having larger cores. 

Our simula t ions largely confirm the results from 
iRead fc Gilmord (|2005|), where a rapid removal of gas re- 
sul ts in a conve r sion f rom cusp to core as stated first 
by iNavarro et al.l (| 19961 ) . As gas cools and flows into the 
halo, the center of the dark matter halo is adiabatically 
compressed. Without star formation, the central gas pres- 
sure builds up, eventually stops further inflow, and even 
makes the gas re-expand somewhat. This re-expansion hap- 
pens rapidly enough for the DM halo to respond non- 
adiabatically: the central DM density experiences a net low- 
ering and the cusp is transformed into a core. With star 
formation turned on, feedback is responsible for a fast re- 
moval of gas from the central parts of the DM halo, with 
the same effect : a conversion from a cu sp to a core. 

Unlike us, iGovernato et al.l (|2010l l found that the den- 
sity threshold for star formation needed to be high enough 
for a cusp-to-core conversion to occur. Only for nsF ^ 
10 cm~^ does supernova feedback lead to sufficient gas mo- 
tions to flatten the cusp in their simulated dwarfs, which 
are taken from a larger cosmological simulation. In contrast, 
in our more idealized, initially spherically symmetric setup, 
even a low density threshold leads to sufficient gas outflow 
for the cusp to flatten. 

3.2 Star formation histories 

In Fig. [21 we show the star- format ion histories (SFHs) of 
different realizations of the NOT mass model. Also, in table 
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— — nsF = 6 CT71 

■ t = 6.00 Gyr 
t = 4.00 Gyr 

■ t = 2.00 Gyr 



10- 



12 3 4 

radius [kpc] 



12 3 4 
radius [kpc] 




Figure 3. The density distribution of the ISM at different times 
for the least massive galaxy, N03, with different density threshold 
and a fixed feedback efficiency of 0.1. 



O the starting time of star formation is tabulated along with 
the final total stellar mass. Several conclusions can be drawn: 

• The delay between the start of the simulation and the 
start of the first star- format ion event is an increasing func- 
tion of nsF- This appears logical: it takes longer for the gas to 
collapse to higher densities and ignite star formation. Com- 
paring different mass models, star formation starts earlier 
in more massive models for a given nsF- This is most likely 
due to the more massive models having steeper gravitational 
potential wells, increasing their ability to compress the in- 
flowing gas. 

• If nsF is increased while epB is kept fixed, more stars 
are formed (e.g. going from the green to the blue curve or 
similarly from the cyan to the magenta curve in Fig. [2]). This 
is because gas collapses to higher densities and the feedback 
is no longer able to sufficiently heat and expel this gas and 
to interrupt star formation. 

• Related to the previous point, the SFR also becomes 
more rapidly varying if nsF is increased while epB is kept 
fixed. The reason is that in the small high-density star- 
forming regions, feedback can only locally interrupt star for- 
mation during short timespans. At lower nsF, star formation 
is more widespread, leading to more global behavior: as su- 
pernovae go off, star formation can be completely halted. 

• Increasing epB while nsF is kept fixed leads to a decrease 
in star formation (e.g. going from the blue to the cyan curve 
in Fig. [2|). This is because once feedback is strong enough, 
it is able to extinguish star formation, even at high gas den- 
sities. 

• The most low-mass models fail to form stars for high 
nsF values. E.g. no stars form in the N03 models for nsF > 
0.1 cm~^. This is due to the masses of these models being 
too small for gas to collapse to densities where stars can be 
formed. This point is further elaborated in the next para- 
graph. 



3.3 Density distribution of the ISM 

In Fig. O the density of the ISM is plotted as a function of 
radius. For the N03 model in the left panel a density thresh- 
old of 0.1 cm~^ was used while for the model in the right 
panel, the density threshold was set to a value of 6 cm""^. 



nsF = 6 cm ^ 
t - 4.00 Gyr 
t = 0.46 Gyr 



5 ^ 




2 3 

radius [kpc] 

Figure 4. The density distribution of the ISM at different times 
for the NOT model, with different density thresholds and a fixed 
feedback efficiency of 0.7. 



The red points show the gas distribution at the moment 
just before the start of star formation in the case of nsF = 
0.1 cm~^. Since up to that moment, all models have expe- 
rienced the same evolution, there is no difference between 
the red points in both panels. As can be seen in the left 
panel, the gas density in this N03 model reaches the star- 
formation threshold and star formation occurs. Moreover, 
the influence of supernova feedback can be seen in the green 
and blue points, where gas expands to larger radii and lower 
densities after having been heated. As is clear from the right 
panel, for nsF = 6 cm""^ the gas simply keeps falling in. It 
will continue to do so during the first 4 Gyr until the built- 
up central pressure causes the gas to re-expand again. No 
stars are formed during the course of this simulation. 

As the density threshold is increased to higher values, 
star formation tends to occur more and more in small col- 
lapsed clumps. This becomes clear when comparing the pan- 
els from Figs. [3] and m The latter shows the gas density dis- 
tributions of two NOT models with nsF = 6 cm~^ and nsF 
= 50 cm""^. While the nsF = 50 cm""^ model only exhibits 
star formation in a small number of discrete high-density 
clumps, the nsF = 6 cm~^ model lacks such well-defined 
clumps and star formation occurs more widespread. 



3.4 Scaling relations 

In this section we discuss the properties of each of our models 
and draw some conclusions regarding the influence of the 
nsF and cfb parameters on the models. An overview of some 
basic properties can be found in Table [2] 



3.4- i Half-light radius Re 

The half-light radius, or effective radius, denoted by Re, 
encloses half of a galaxy's luminosity. In panel a.) of Fig. \5\ 
Re is plotted as a function of the V-band magnitude. The 
following trends can be observed in this figure: 

• For a fixed nsF, the effective radius varies only very 
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Table 2. Final properties of our large set of simulations. Columns: (1) model number (see Table[T]), (2) density treshold for star formation, 
(3) feedback efficiency, (4) final stellar mass, (5) starting time of star formation, (6) half-light radius, (7) mean surface brightness within 
the half-light radius, (8) central one dimensional velocity dispersion, (9) mass-weighted metallicity, (10) central surface brightness, (11) 
Sersic parameter, (12) circular velocity. 



Model 


nsF 


epB 


M^,f 


ATsF 


Re 


/e 


CriD,c 


V -I 


[Fe/H] 


MO 


n 


Vc 




[cm-3] 




[106 Mo] 


[Gyr] 


[kpc] 




[km/s] 






[mag] 




[km/s] 


N03 


0.1 


0.1 


0.285 


0.342 


0.100 


3.1T6 


6.806 


0.835 


-1.183 


23.603 


1.496 


16.1T6 


N05 


0.1 


0.1 


5.66T 


0.168 


0.230 


10.TT9 


12.190 


0.860 


-1.236 


23.1T0 


0.959 


20.396 


N05 


6 


0.1 


1T.86T 


0.546 


0.130 


88.566 


12.253 


0.910 


-0.659 


20.920 


1.025 


22.5T4 


N05 


6 


0.3 


4.049 


0.546 


0.142 


18.222 


9.131 


0.8T0 


-1.130 


23.005 


0.8TT 


20.483 


N05 


6 


0.5 


2.021 


0.546 


0.134 


11.813 


8.310 


0.839 


-1.302 


22.586 
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slightly throughout the epB -range and this without a clear 
trend between Re and epB- However, for a fixed nsF and 
dark-matter mass the stellar mass and consequently the lu- 
minosity decrease with increasing cfb- This is due to star 
formation being shut down more rapidly when feedback is 
more effective. As a result, galaxies tend to have higher stel- 
lar densities for smaller cfb- 

• For a fixed epB of 0.1, an increase of nsF from 0.1 to 
6 cm""^ results in a decrease of the effective radius. This is 
due to the size of the region where the SFC are fulfilled, 
which is much smaller for nsp = 6 cm""^ than for nsF = 
0.1 cm~^, and the feedback is too weak to overcome this. In 
the case of an increase of nsF from 6 to 50 cm""^, the effective 
radius increases which is caused by the higher star forma- 
tion peaks resulting in more supernovae explosions which 
redistribute the gas more efficiently. 

• The simulations with high density threshold, nsF > 
0.1 cm~^, and high feedback efficiency, cfb > 0.1, have ef- 
fective radii which are in agreement with the observations. 

From this scaling relation we can constrain the cfb- 
parameter to he higher then 0.1 to produce galaxies with ef- 
fective radii in agreement with observations of dwarf galax- 



3.4-2 The fundamental plane 

The fundamental plane (FP) is an observed relation between 
the effective radius, i?e, the mean surface brightness within 
the effective radius, /e, and the central velocity dispersion, 
Gc of giant elliptical galaxies. It is a linear relation, given by 

log(i?e) = -0.629 - 0.845 log(/e) + 1.38 log(crc), (10) 

between the logarithms of these quantities (|Burstein et al.l 
ll997l V In panel b.) of Fig.O we plot the "vertical" deviation 
of the simulated galaxies from the giant galaxies' FP. 

Dwarf galaxies generally lie above the FP in this pro- 
jection. This is thought to be a consequence of their having 
shallower gravitational potential wells than giant galaxies. 
This, together with the feedback, results in more diffuse sys- 
tems. Models with a high star- format ion threshold in com- 
bination with a low supernova feedback turn out to be very 
compact. They actually populate the FP at low luminosities. 
However, this region of the three-dimensional space spanned 
by log(i?e), log(/e), and log(crc) is observed to be devoid of 
galaxies. Hence, models with low stellar feedback, cfb up to 
0.3, and high density thresholds, nsF > 0.1 cm~^, can be 
rejected. 

3.4.3 ColorV-I 

Fig. O panel c.) shows the V — I color in function of the 
y-band magnitude. The color scatter between the differ- 
ent models is rather small. The observed galaxies follow 
a mass-metallicity relation so the metallicity generally in- 
creases with the galaxy (stellar) mass, resulting in increas- 
ing V — I values for increased galaxy mass. Within the rela- 
tively small mass range covered by the models, color is only 
a very weak function of stellar mass. For a fixed feedback ef- 
ficiency, when increasing the density threshold the V — I also 
increases slightly resulting in bluer galaxies for the models 
with low density threshold. This is due to the effect that 



stars are formed in more metal enriched regions in the mod- 
els with high density threshold. When the density threshold 
is kept constant and only the feedback efficiency is increased 
the V — I slightly decreases, so the models get slightly bluer 
due to a dilution of the gas when it is more spread out by 
supernovae explosions. 



3.4.4 Metallicity 

In panel d.) of Fig. [5] a plot of iron content [Fe/H] as a func- 
tion of the V-band magnitude is shown. The mass- weighted 
value of [Fe/H] is a measure of the metallicity of a galaxy. 
The yellow and magenta dots represent observational data 
from dwarf spheroidal and dwarf elliptical galaxies and irreg- 
ular dwarf galaxies, respectively. Some general conclusions 
we can take away from this figure are: 

• Low-mass models with low density threshold, nsF ~ 
0.1 cm""^, and low feedback, cfb ~ 0.1, keep forming stars 
throughout cosmic history and do not expel enriched gas. 
As a consequence, they turn out to be too metal rich, com- 
pared with observed dwarf galaxies. Models with higher nsF 
compare much more favorably with the data in this respect. 

• For a fixed nsF, increasing epB, produces more metal 
poor galaxies. This is likely due to the fact that the in- 
creased feedback extinguishes star formation more rapidly 
and disperses the metal enriched gas more widely. 

• Increasing nsp at fixed cfb and fixed mass, results in 
an increase of the metallicity and of the stellar mass when 
going from nsF = 0.1 cm~^ to nsF = 6 cm~^. A further 
increase of nsF at fixed cfb, up to nsF = 50 cm""^, has a 
much smaller impact on metallicity and stellar mass. The 
former is likely due to more vigorous star formation in less 
easily dispersable high density regions. 



3.4-5 Surface brightness profiles 
We fitted a Sersic profile, of the form 



I{R) = Joe" 



^^y/- 
.^0; 



(11) 



to the surface brightness profiles of the simulated galaxies. 
The Sersic parameter n and the central surface brightness 
/io are plotted respectively in the panels e.) and f.) of Fig. [5] 
as a function of the l/-band magnitude. 

• For a fixed nsF , when increasing the cfb , there is a weak 
trend for the Sersic parameter n and the central surface 
brightness to decrease. More vigorous feedback appears to 
result in more diffuse dwarf galaxies, as one would expect. 

• As an echo of the Re — My relation, simulations with 
high density threshold, nsF > 0.1 cm~^, and low feedback 
efficiency, cfb = 0.1 — 0.3, are systematically too compact, 
with /io ^ 20 mag arcsec"^, compared with the observations. 

• The models with high density thresholds and strong 
feedback are in general agreement with the observations. 



3.4-6 The Tully-Fisher relation 

Panel a.) of Fig. [6] shows the B-band Tully-Fisher relation 
(TFR) between the circular velocity, denoted by 14, and the 
luminosity in the B-band, Lb- The simulations are compared 
with observational data and with the Tully-Fisher relation 
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Figure 5. Some scaling relations and the surface brightness parameters as a function of the magnitude. In a.), the half-light radius Kq is 
plotted, b.) shows the vertical deviation of the simulated dwarf galaxies from the giant galaxies' FP, in c.) the V — I color is plotted, d.) 
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Figure 6. The top panel shows the Tully-Fisher relation between 
the circular velocity and the luminosity in the B-band. The full 
gray line shows the TF relation for early type galaxies, the dashed 
gray line is the TF re lation of spiral galaxies as determined by 
IPe Riicke et al.l (|2007l V The lower panel shows the Faber- Jackson 
relation between the velocity dispersion and the luminosity in the 
B-band. 



for early-type (full gray line) and fo r spiral galaxies (dotted 
gray line) that was determined by IPe Riicke et al.l (|2007l ) . 
All simulations predict that the TFR becomes substantially 
shallower in the dwarf regime, below luminosities of the or- 
der of Lb ^ 10^ Lq^b- This can be seen as a consequence 
of the very steep Mstar — Mhaio relation in the dwarf galaxy 
regime (see paragraph I3.4.8P . For a fixed nsF , an increase 
in feedback efficiency does not influence Vc very much since 
there are so few stars that Vc is set by the dark-matter halo. 
The effect on the stellar mass, and consequently on Lb, is, 
however, quite large. Therefore, increasing cfb at fixed usf 
and dark-matter mass causes galaxies to shift leftwards in 
panel a.) of Fig. [S] Except for this effect, once nsF and cfb 
are raised above their minimum values of 0.1 cm""^ and 0.1, 



respectively, there is no significant differences between the 
TFRs traced by the different series of models. 



3.4-'^ The Faber- Jackson relation 

The Faber- Jackson (FJR) relation, plotted in panel b.) of 
Fig. [6] is the relation between the stellar central velocity dis- 
persion and the luminosity in the B-band. The stellar central 
velocity dispersion is a projection of the velocity dispersion 
along the line of sight. This is measured by fitting an expo- 
nential function to the dispersion profile and retaining the 
maximum of the function as the central value. 
From this figure we see: 

• For a fixed nsF, when increasing the cfb, the velocity 
dispersion decreases first after which it settles around a value 
which depends on the dark-matter mass of the model. 

• For a fixed cfb, when increasing nsF, only a minor in- 
fluence on the velocity dispersion is observed. 



3.4 'S The Mstar -Mhaio relation 

In Fig. [71 the Mgtar-Mhaio relation of the simulations at 2: = 
is plotted. We can make similar conclusions here as were 
made in the SFH section: 

• If nsF is fixed, the stellar mass will decrease if the cfb 
is increased. This is what was expected because with more 
feedback the gas is distributed over a larger area and the 
infall of the gas to the appropriate density threshold will 
take longer. 

• If Cfb is fixed, for increasing nsF, the stellar mass in- 
creases too. When feedback is very small, the gas density will 
stay high and the star formation will not be interrupted, re- 
sulting in a high stellar mass. The effect is smaller for higher 
feedback. 

In Fig. [71 our different sets of models are found to be in 
agreement with the results from the Aquila simulation where 
a density threshold of 10 cm""^ and a feedback efficiency of 
0.7 was used. While the initial conditions of our dwarf galaxy 
simulations are admittedly quite simplified, they do have 
high spatial resolution and realistic implemented physics. It 
is therefore encouraging that they compare favorably with 
cosmological simulations like the Aquila simulation, which 
have cosmologically well motivated initial conditions but in 
which dwarf galaxi es are very close to the resolution limit 
(jSawala et al.ll201ll ) . However it is impossible by further tun- 
ing of the feedback efficiency and/or the density thre s hold t o 
reproduce the trend that was derived bv JGuo et al.l (|2010h . 

By increasing the density threshold and feedback ef- 
ficiency, the stellar mass is reduced by almost two orders 
of magnitude, but there still remains a difference of many 
orders of magnitude between our sim ulations and the M^- 
Mhaio relation from IGuo et al.l (|2010|). It is also interesting 
to notice that although our models do not reproduce the 
relation, they do have a very similar slope. 



© 2011 RAS, MNRAS 0Q0.fTIfT2l 



10 A. Cloet-Osselaer et al. 



1012 

10" 

lOio 

I 10« 

I 108 

10^ 
10*5 



10^ I' 
10^ 



Sawala et al. 2011, isolated galaxies 
Sawala et al. 2011, satellite galaxies 
Mandelbaum et al. 2006 
nsF = 0.1 cm^"^, epB = 0.1 
nsF = 6 cm~^ 
nsF = 50 cm^"^ 




0.9 
0.7 
0.5 
0.3 
0.1 



0.9 

0.7 

( 

0.5^ 
0.3 



10^ 



lO^o lO^i 

Mdm [Mo] 



10 



12 



10 



13 



10^ 



Figure 7. The stellar mass versus the DM halo ma ss , plot - 
ted in comparison with the models by I Sawala et al.l (|201lh . 
The gray dots show d ata from gravitational lensing from 
iMandelbaum et al.l (|2006r ) . The black line is the trend for this 
relation that was determined bv lGuo et al.l (|2010l ). 



4 DISCUSSION AND CONCLUSIONS 

4.1 Cusp to core 

Whether the halo density profile is cusped or cored has been 
a point of discussion for quite some ti me. Observationally , 
evidence for cored DM profiles is found (| Gentile et al .11 20041 ) , 
but from cosmol ogical DM sim ulations a cusped densit y pro- 
file is deduced (jNa varro etHI FlOQe: Mo ore et al.lll996l ). The 
inherent limitation due to the angular resolution of the ob- 
servations i s ruled as a cause of th e observed fiat d ensity 
profiles by Ide Blok k Bosmal (l2QQ2l ). iGentile et all (|2005l ) 
also excluded the possibility of non-circular gas motions 
which might result in a rotation curve that is best fitted 
by a cored halo, while the dark matter halo actually has a 
cuspy profile. However, from the simulation point of view, 
iMashchenko et al.l (|2006l l mentioned a natural transition of 
a cusp to a flattened core when the dark matter halo is grav- 
itationally heated by bulk gas motions. 

Our simulations are set up with a cusped NFW halo 
in agreement with cosmological simulations. The infall of 
gas causes an adiabatic compression of the dark halo. When 
gas is evacuated from the central regions, be it by a fast 
re-expansion as the gas pressure builds up or by supernova 
feedback, the dark-matter halo reacts non-adiabatically and 
kinetic energy of the gas is transferred to the dark matter. 
This results in a flattening of the central density and so 
the cusp is converted into a core. We can conclude that 
the conversion of the cusped halo density profile to a cored 
profile is realized by the remov al of baryons from the galaxy 
center (iRead fc Gilmord 12005"), whether this is due to a re- 
expansion of the gas or by feedback effects or by another 
process. 



4.2 Degeneracy 

By increasing both the density threshold and the feedback 
efficiency, the simulated galaxies move along the observed 
kinematic and photometric scaling relations. These two pa- 
rameters, the feedback efficiency cfb and the density thresh- 
old nsF , correlate with each other and an increase of the one 
can be counteracted by an increase of the other, resulting 
in galaxies with similar properties. To be more specific: the 
individual galaxies are drastically different for different pa- 
rameter values but they all line up along the same scaling 
relations and can therefore be seen as good analogs of ob- 
served dwarf galaxies. 

The feedback efficiency quantifies the fraction of the 
10^"^ ergs of energy that are released during a supernova ex- 
plosion and thermally injected into the ISM. For each value 
of the density threshold we can determine the feedback effi- 
ciency range for which the models are in agreement with the 
observations, although we are not able to deduce a unique 
^sr/eFB-combination which would be the "correct" repre- 
sentation of the physical processes that happen in galaxies. 

For a certain density threshold, a lower limit of the cor- 
responding epB -parameter can be determined from the effec- 
tive radius: the galaxies become too centrally concentrated 
when the feedback is too low. From the scaling relations we 
cannot deduce an upper limit for the cfb -parameter, but 
one could argue that the ISM cannot receive more energy 
than there is released by the supernova explosion, resulting 
in a maximal value for the feedback efficiency of 1. 

In the case of a density threshold of nsF = 0.1 cm""^, 
the models are generally in good agreement with the ob- 
servations besides the somewhat high metallicities. This is 
also the reason why the feedback efficiency was not varied 
in this case. If we compare the high density threshold mod- 
els, nsF > 0.1 cm~^, with the observations we can conclude 
that the feedback efficiency should be larger then '^ 0.3. For 
a density threshold of nsF = 6 cm~^, we prefer a value of 0.7 
for the feedback. Similarly we prefer a feedback efficiency of 
0.9 in the case of a density threshold of nsF = 50 cm""^ 

The fact that different tt-sf/cfb -combinations result in 
simulated galaxies with properties that are in agreement 
with the observations invokes a warning for future simu- 
lations and indicates that there is still some work left to 
determine the density of the star forming regions and the 
fraction of supernova energy that is absorbed by the ISM, 
quantities which are hard to determine observationally. 

There are however other parameters that might influ- 
ence the starformation rate and our degeneracy, which are 
not investigated here: 

• Given the fact that the star- format ion efficiency c* was 
found by other authors not to have a significant impact on 
stellar mass, we did not investigate it in detail in this paper. 

• The choice of the IMF, for which in our simulations a 
Salpeter IMF is used, determines the mass distribution of 
stars. The fraction of high-mass stars influences the number 
of SNIa and SNII explosions and as a consequence it will 
influence the amount of feedback and the chemical evolution. 
However, given the large number of IMF parameterizations 
available in the literature, testing them is a very daunting 
task which falls outside the scope of this paper. Moreover, 
part of the IMF- variation is quantified approximately by the 
variation in cfb which we do investigate. 
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density threshold com pared to the relati on of IGuo et ahl l2010l . 
other simulations from [ Sawala et al.ll201ll and observations from 
iMandelbaum et al.llSooa 



• There are other possible feedback implementations, 
next to the release of feedback energy as thermal energy to 
the gas. It also could be released as kinetic en erg:y by kick- 
ing t he gas particles or by blast- wave feedback ([Maver et alJ 
I2OO8I I. 

• Other implementations of star formation, e.g. based on 
a subgrid model of H2-formation (;Pelupessv et al.M2QQ4 ). are 
possible. 



- Not properly taking into account the effects of reion- 
isation may lead to an overestimation of the gas content 
of dwarfs and an underestimation of the gas cooling time. 
However, even taking into account reionisa tion, the dwarf 
galaxies simulated by I Sawala et al.l (|201lh had much too 
high stellar masses. 

- At a given gas density, the star- format ion efficiency 
of dwarf galaxies could be lower than that of more mas- 
sive stellar systems because of their lower metallicity and 
hence lower dust content. This could be mimicked by re- 
ducing the star- format ion efficiency parameter c^ ( see ea . 
dS])) in the dwarf regime. However, IStinson et al.l (|2006|) 
have shown that, because of self- regulation, the star- 
formation rate is very insensitive to this parameter: vary- 
ing c^ between 0.05 and 1 left the mean star- format ion 
rate virtually unchanged. 

- External processes such as ram-pressure stripping 
and tidal stirring may lead to a premat ure cessation 
of st ar formation and hence lower Mstar rM aver et al.l 
1 20061 ). However, these processes are only effective if the 
gravitational potential wells of dwarf galaxies are suf- 
ficiently shallow and if they are stripped early enough 
in cosmic history, before they converted their gas into 
stars. It is uncle ar whether these constraints are met. In 
iDe Riicke et al.l ([2010 ), and references therein, it was ar- 
gued that the number of red-sequence, quenched dwarf 
galaxies increased significantly over the last half of the 
Hubble time and that the dwarf galaxies now residing in 
the Fornax cluster were accreted less than a few crossing 
times age (i.e. less than a few Gyr). This timescale would 
have left dwarf galaxies ample time to form stars before 
entering the cluster. 



4.3 The dwarf galaxy dark- matter halo occupancy 

To conclude. Fig. [8] shows the models which best agree with 
the observations for each density threshold that was used 
in our analysis. Increasing nsF together with cfb leads to a 
strong reduction, of almost two orders of magnitude, of the 
stellar mass, especially in the most massive models. How- 
ever, with the physics included in our simulatio ns, we are 
unabl e to reproduce the Mstar — Mhaio relation of lGuo et al.l 
(|2010l V Surprisingly, the best models trace a Mstar — Mhaio 
relatio n with a s l ope t hat is similar to that of the rela- 
tion of IGuo et al.l (|2010l V Our simulations are in agreement 
with results from cosmological simulations, which have, how- 
ever, much lower spatial resolution in the dwarf regime 
I Sawala et al.l (| 20111 ) . We did not explore yet higher values for 
nsF and epB because it is clear from Fig.[8]that the reduction 
of Mstar stagnates for high nsF-values. Moreover, to compen- 
sate for the high density threshold, an unphysical large value 
for epB, higher than 1, would be required. Thus, we arrive at 
(nsF = 6cm~^, epB ^ 0.7) and (nsF = 50cm~^, epB ^ 0.9) 
as the models which are in best agreement with the observed 
photometric and kinematical scaling relations and with the 
Mstar — Mhaio relation derived directly from cosmological 
simulations. 

While it appears impossible to place isolated dwarf ga lax- 
ies on the Mstar — Mhaio relation of IGuo et ahl (|2010l l , it 
is possible to envisage external influences that may further 
reduce Mstar, as already mentioned in the Introduction: 
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